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CHANNEL ESTIMATION FOR OFDM SYSTEMS 



The present invention relates generally to methods of channel estimation in 
wireless Orthogonal Frequency Division Multiplexing (OFDM) systems, and in 
5 particular to methods of channel estimation using Linear Minimum Means Square 
Error (LMMSE) estimation techniques. 

Orthogonal Frequency Division Multiplexing (OFDM) is a high spectral 
efi&ciency type of multi-carrier modulation system, which has many advantages of 
single carrier systems, especially for high data rate transmission in time dispersive 
10 channels. Transmitted diversity is an effective method to further improve vrareless 
conmiunication systems in fading environments. Space-time coded OFDM systems 
with transmitter diversity capable of reliable high data rate wireless communications 
promise to be an effective alternative for broadband wireless services. However, 
space-time coded systems require accurate estimation of channel frequency responses. 

IS Traditional one-dimensional channel estimation techniques for OFDM systems 

include (a) Leased Squares (LS), (b) Minimum Means Square Error (N4MSE) and (c) 
Linear Minimum Means Squared Error (LMMSE) estimation techniques. LS 
estimators have low complexity, but suff^ from a high Means Square Error (MSE), 
especially if the system operates with low signal to noise ratios. On the other hand, 

20 MMSE estimators, based on time-domain channel statistics, are hig^y complex and 
require significant numbers of multipli^ and adders in any practical implementation. 
MMSE estimators provide good performance for sample spaced chaxmel 
environments, but have limited performance for non-sample spaced channels and higfh 
signal to noise ratios. 

25 LMMSE estimators provide good performance for sample spaced and non- 

sample spaced chaonels. Nevertheless, practical implementations of LMMSE 
estimators suffer from being highly complex and require a large number of 
computations to be performed in order to achieve accurate channel estimatioiL 



30 an OFDM system with transmitter diversity that is simple and efficient, and minimises 
the computational complexity of existing channel estimation techniques. 



It would be desirable to provide a method for performing channel estimation in 
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It would also be desirable to provide a method for performing chaiuiel 
estimation that alleviates or overcomes one or more problems of known channel 
estimation techniques. 

One aspect of the present invention provides a method for performing linear 
5 channel estimation in an orthogonal frequency-division multiplexing systeni, the 
method including the steps of: 

receiving transmitted pilot symbols from a plurality of transmit antennas; 

fomiing a least-squares estimation matrix from the transmitted pilot symbols; 

foradng a sparse smoothing matrix approximating a fixed weighting matrix, 
10 wherein each row vector in the sparse smoothing matrix contains one or more of the 
strongest weights in each row of the fixed weighting matrix; and 

deriving a channel estimation matrix Sxmi the sparse smoothing matrix and the 
least-squares estimation matrix. 

In one embodiment the sparse smoothing matrix is defined according to: 

where Ej{k) is the row energy of the sparse smoothing matrix with non-zero terms 
Wj{k,m) formed from the M strongest weights of the k*th row of the fixed weightiQg 
matrix Wj{k\ k represents the frequ^cy bin numba- and j the transmitting antenna 
number. 



15 



Wj{k) 



20 The repeated pilot symbols may be preceded and/or followed by a cycKc prefix 

and may be transmitted on interleaved sub-carriers from the plurality of transmit 
antennas. 

Alternatively, flbie independent pilot symbols, may each be preceded and /or 
followed by a cyclic prefix, and may be transmitted on interleaved sub-carriers from 
25 the plurality of transmit antmnas. 



la another alternative, each pilot symbol may be preceded and/or followed by a 
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cyclic prefix that is transmitted on interleaved sub-carriers fiom the plurality of 
transmit antennas. 

Preferably, a cyclic prefix window length, or delay spread approximation length 
is chosen to enable &e real and imaginary parts of the fixed weighting matrix to 
5 contain equal or zero entries. The lengfli of the cyclic prefix window or the delay 
spread approximation can be (l+N/2) or (l+N/4) where N is the length of the Inverse 
Discrete Fourier Transform used to form the pilot symbol. 

In a preferred arrangement the step of forming a sparse smoothing matrix 
includes: 

1 0 calculating a plurality of possible sparse smoothing matrices; 

storing the plurality of matrices in a storage device; and 

selectively retrieving one of the plurality of possible sparse smoothing matrices 
firom the storage device. 

The storage device may conveniently be a look-up table. 

15 The smoothing matrix may be selected for retrieval firom the storage device 

according to characteristics derived firom the least squares estimation matrix. 

The characteristics may include any one or more of the signal to noise ratio 
SNR, the root mean square delay spread of the power delay profile r„ and the delay 
spread of the power delay profile r x- 

20 The method may fiirth^ include the step of: 

making coefBcients of ttie fixed weighting matrix real by performing a cyclic 
shift to locate the channel impulse response symmetrically around zero. 

Conveniently, cyclic shift may be performed in either the time domain or by an 
equivalent linear phase rotation in the fi-equency domain. 

25 The method may fiirther include the step of: 

using a symmetrically shaped delay spread approximation for the channel 
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estimatioii. The delay spread approximation may be rectangular-shaped. 

Anotiier aspect of the invention provides a channel estimator for use in an orthogonal 
frequency-division multiplexing system, the channel estimator including: 

a least-squares estimation unit for forming a least-squares estimation matrix 
5 from pilot symbols transmitted from a plurality of transit antennas; 

a matrix foraiation unit for forming a sparse smoothing matrix approximating a 
fixed weighting matrix, wherein each row vector in the sparse smoothing matrix 
contains one or more of the strongest weights in each row of the fixed weighting 
matrix; and 

10 a channel estimation unit for forming a channel estimation matrix from the 

sparse smoothing matrix and the least-squares estimation matrix. 

Convenientiy, the matrix formation unit may include: 

a storage device for storing a plurality of possible sparse smoothing matrices; 

and 

15 a matrix selection unit for selectively retrieving one of the plurality of possible 

sparse smoothing matrices from the storage device. 

The storage device may be a look-up table. 

The matrix formation unit may act to select the sparse smoothing matrices for 
retrieval from the storage device according to characteristics derived from the least 
20 • squares estimation matrix. 

In order to assist in arriving at an understanding of the present invention, a 
preferred embodiment is illustrated in the attached , drawings. However, it should be 
understood that the following description is illustrative only and should not be taken in 
any way as a restriction on the generality of the invention as described here above. 

25 In the drawings: 

Figure 1 is schematic diagram of an OFDM system; 
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Figure 2 is a schematic diagram of a chaxmel estimator forming part of a 
receiver in the OFDM system of Figure 1 ; 

Figure 3 is a flow chart illustrating operation of the channel estimation of 
Figure 2; 

5 Figure 4 is a diagrammatic representation of three diflferent pilot symbol 

allocation schemes for use in the channel estimation process shown in Figure 3; 

Figure 5 is a diagrammatic representation of the symmetrical location around 
zero of the channel impulse response and uniform delay spread approximation used in 
the LMMSE channel estimation shown in Figure 3; 

10 Figure 6 shows the mean squared error peiformance vs complexity of the SWC 

method compared to the SVD method; and 

Figure 7 shows the mean squared error performance vs SNR for the SVD and 
SWC schemes. 

Referring now to Figure 1, there is shown generally an OFDM based system 10 
15 which exploits channel estimation and signal detection operations in equalisation. A 
digital signal source 12 is protected by channel coding jfrom a diannel encoder 14 and 
is interleaved by an interleaver 16 against fading phenomenon. After this, the binary 
signal is modulated by an OFDM modulator 18 and transmitted over a multipatfa 
fading channel 20. During transmission, noise 22 is added. 

20 The simi signal is received at a receiver filter 24, which can take the form of a 

DFT (Discrete Fourier Transform), and the output of the filter then passed to a signal 
detector 26, Due to the multipath channel transmission, some inter-symbol 
interference occurs in the received signal. Accordingly, the signal detector 26 requires 
knowledge of the Channel Impulse Response (CIR) characteristics in order to ensure 

25 successfiil removal of the inter-symbbl interference. The channel impulse response 
characteristics are determined by a channel estimator 28. After detection, the signal is 
de-iaterleaved by a de-intCTleaver 30 and the chaimel decoded by a channel decoder 32 
to extract the original message. 

Transmitter diversity is achieved in the OFDM system 10 shown in Figure 1 by 
30 flie use of multiple transmit antennas. To enable channel estimation, pilot symbols are 
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simultaneously sent from tbe multiple transmitter antennas on interleaved sub-carriers. 
At the receiver end, the LMMSE channel estimator 28 identifies channel 
characteristics in the non-measured sub channels by interpolating different sets of 
measured sub channels from each specified antenna. 

5 In a downlink diversity environment with two transmit antennas and one 

receiver, the two transmit antennas j = 1, 2 simultaneously send to OFDM pilot 
symbols on K interleaved sub-carriers. Tbe pilot symbols X, and X2 are defined as 
follows: 

xl = {ao, 0, ai, 0, aa, aK/2-1, 0} 

10 x2 = {0, bo, 0, b,, 0, b2, 0, bK/2-1} (1) 

where a* and bk are arbitrary complex numbers with magnitude of 1 . 

Each of these signal forms an OFDM block. With the channel impulse 
response confined to a cyclic prefix (CP) length, the Digital Fourier Transform (DFT) 
of the received symbols can be given by 

15 y(k)^j;^Hjik)xjik)+vik) (2) 

where k = 0, 1, K - 1 denotes the sub-carrier number, Hj(k) is the channel 
frequency response corresponding to transmit antenna J and v(k) is the additive 
complex Gaussian noise with zero mean and variance one. 

In this exemplary embodiment, the channel estimator 28 is a packet-type 
20 channel estimator, where only the firequency correlation of the channel is used in the 
channel estimation. The firequency domain correlation depends on the multipath 
channel delay spread and can be described by a frequency domain correlation fimction 
r f (k). For an exponentially decaying multipath power delay profile, the frequency 
domain correlation fimction rf (k) can be given by 

25 rf(k) = ^ — 7—^ (3) 

where is the root-mean square (rms) delay spread of the power delay profile and 
4r denotes the sub-carrier spacing. 
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The LMMSE channel estimation vector corresponding to the /h 
transmitter in a 2 x 1 diversity system can be obtained as follows: 

^y=W?>/y (4) 

where R„ ^ —Rmp and iJs s = | Rpj> + ] are the correlation matrices of size Kx 

5 K/2 and KJ2 x KJ2 respectively [3]. I is the identity matrix and SNR is fhe expected 
value of SNR. Pj is the least-squares (LS) estimation vector of length K/2 at the pilot 
positions corresponding to antenna j, given by 

Pj^x^yj (5) 

where Xj is a diagonal matrix containing the transmitted pilot points xj{k) given by 
10 (1). 

The best low-rank approximation of R^h^J^^^J^j^ given by Singular Value 

Decomposition (SVD). Then, with the appropriate substitutions in (4), the rank- r 
estimator is defined by 

15 where and Vf are unitary matrices, and is the r x r iq^per left comer diagonal 

matrix, containing the strongest singular values. The siq^erscripts (.)'' and (.) ^ denote 
rank-r and Hennitian transpose respectively. 

In channels wi& large delay spreads, the rank-r approaches a value ofK/S^ the 
low rank approximation no longer reduces the estimator complexity. 

20 The chamiel estimator 28 provides an alternative sparse approximation of the 

fixed weighting matrix, namely LMMSE by significant weight catching (SWC). For 

notional convenience, the equation (4) can be rewritten. 

• . 
Hj'^WjPj (7) 
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where Wj = R^^^^RZL is the fixed weighting matrix (otherwise known as the 

interpolation matrix). 

Some row entries of the Wj contain stronger weights than the others, with the 
strongest values on its diagonal. 

5 The channel estimator 28 acts to restrict the frequency domain of titie fixed 

weigjiting matrix Wj to be a sparse (i.e only including limited number of non-zone 

elements) smoothing matrix containing the M strongest weights in each row, where 
M The sparse smoothing matrix approximating the fixed weighting matrix is 

obtained fix>m: 

10 Ej{k) = argm^ll^^^ (8) 

where Wj{k) denotes a row vector firom the fixed weighting matrix. 

Figure 2 shows a practical implementation of the channel estimator 28. A de- 
multiplexer block 40 acts de-interleave pilot symbols into streams based on the 
transmit antenna firom which the pilot symbols originated. Least squared estimators 42 

15 and 44 are based on known pilot data and receive the pilot symbol streams fi-om the 
de-multiplexer block 40. Inverse Fast Fourier Transform (IFFT) blocks 46 and 48 act 
to estimate the impulse response fi*om which the route mean square delay spread (in 
blocks 50 and 52) and signal-to-noise ratio estimates (in blocks 54 and 56), together 
with other features, for example the absolute delay spread, are extracted. A common 

20 logic block 48 receives the sigual-to-noise ratio estimates and route mean squared 
delay spread estimates and other features, and acts to select an appropriate sparse 
smoothing matrix firom a lookup table stored in the non-volatile memory device 60. 
Rotators 62 and 64 act to rotate the least squared estimates generated by blocks 42 and 
44, which are then multiplied and summed with the sparse smoothing matrix identified 

25 by the common logic block 58, by means of the multiply and sum blocks 66 and 68. 
The rotator block 62 and 64 perform a channel impulse response rotation in the 
fi-equency domain. The multiply and sum blocks 66 and 68 act to smooth and 
interpolate the least squared estimates exploiting the significant weight catching 
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technique of the present inventiorL The rotating blocks 70 and 72 then act to de-rotate 
the output of the multiply and sum blocks 66 and 68 in order to generate the chaimel 
estimates. It should be noted that the de-rotation blocks 70 and 72 can be avoided if 
the data is pre-rotated. 

5 The steps carried out by Ihe chaxmel estimator are depicted in Figure 3. This 

figure shows that initially, at stqp 80, transmitted pilot symbols are received ficom the 
multiple transmit antennas used in the OFDM system with transmitted diversity shown 
in Figure 1. At step 82, the least squares estimation matrix ^is conoputed by the 

channel estimator 28 according to the expression Pj = Xj^yj . 

10 Tlie LMMSE channel estimation efTector Hj can be obtained jfrom the product 

of a sparse smoothing matrix and the least squares estimation. In order to further 
nrnmnise channel estimator coiiq)lexity and improve the estimation accuracy of the 
channel estimator 28, a number of possible sparse smoothing matrices may be 
calculated and stored in a lookup table within the channel estimator 28 beforehand. 

15 In order for this to occur, a channel impulse response is initially obtained by 

performing an hiverse Fast Fourier Transform (IFFT) operation at step 84 on the least 
squares estimation matrix. From the hiverse Fast Fourier Transform, the signal to 
noise ratio, the mean square delay spread of the power delay profile and delay spread 
of the received pilot symbols are firstly calculated. The power delay profile is the 

20 output of the IFFT and it is confined to the length of the cyclic prefix. A noise estimate 
can be taken fi'om the other outpirts to form an SNR estimate. The time between the 
first and last significant multipath component of the power delay profile is the delay 
spread and the rms delay spread can be obtained from: 




25 where the is the amplitude and r,. is the delay of the i'th midtipath component. 

With the knowledge of the aforementioned channel impulse response 
characteristics having been estimated at step 86, the most appropriate interpolation or 
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sparse smoothing matrices is then selected by the channel estimator 28 from a lookup 
table, at step 88. 

At step 90, the LMMSE channel estimation is carried out by computing the 
product of the sparse smoothing matrix selected by the channel estimator at step 58 
5 and &e least squares estimation matrix as determined in step 82. Broadband Wireless 
Local Area Networks (WLANs) incorporate two long OFDM pilot symbols at the 
beginning of a data packet, to enable chaimdi estimation. The pilot symbols are 
preceded by a double length Cyclic Prefix (CP) to effectively eliminate inter-symbol 
interference and inter-carrier interference due to a fading chamieL The foUowing 

10 modified pilot schemes that enable the inclusion of transmitter diversity or multiple 
input multiple output systems within existing OFDM standards have been found to be 
particularly suitable for use with the present invention. The first scheme, shown in 
Figure 4(a) consists of a standard pilot system in which two repeated (in this case long) 
pilot symbols 190 and 102 are preceded with a cyclic prefix 104. In tibis case, the 

IS cyclic prefix is a double length cyclic prefix of 1600 115. 

The second scheme, shown in Figure 4(b) splits titie two repeated pilot symbols 
into two independent pilot symbols 106 and 108» each of which is preceded with a 
cyclic prefix, in this case a single cyclic prefix of length 800 ns. The cyclic prefix 
preceding the pilot symbol 106 is referenced 1 10 in Figure 4, whilst the cyclic prefix 
20 preceding Che pilot symbol 1 08 is referenced 1 12. 

The third scheme shown in Figure 4(c) transmits a single pilot symbol 114 
preceded by a cyclic prefix, in this case a double length cyclic prefix length of 1600 its 
referenced 116 over twice the number of sub channels but half the bandwidth of the 
two previously mentioned sch^es. 

25 The three exemplary schemes ^own in Figure 4 are 4x1 antenna diversity 

system. The first two schemes form two consecutive OFDM pilot symbols Xj (i), i = 
(0, 1) for each antenna J = (1, 2, 4). The third scheme forms only one pilot symbol 
xj (0, i = 0 for each antenna J. All tiiree schemes have a preamble length of 8 jiis. 
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In channels with a limited mobility, the least squares estimation matrix Pj of 

the two repetitive OFDM symbols in the first pilot scheme, shown in Figure 4(a) can 
be obtaxaed in step 82 as follows: 

. 5 where Xj = Xj (i), i = (0,1) is a diagonal matrix of size K/Q x K/Q containing the 
transmitted pilot points xj (&). 

Ihe least Pj squares estimation matrix in the second pilot scheme, shown in 
Figure 4(b) , can be obtained in step 62 by: 

P^=:P^(0)U^(1) (10) 

10 where P, (z) is the LS estimates vector of length K/Q^ coirespondiug to the zth 
received pilot OFDM symbol firom transmitter y, given by: 

p,(f)=Ar;»(/)y,(i) (11) 

Equation (11) also represents the LS estimation vector = ij(/)> i ^ 0 of 

length 1K/Q for the tiiird pilot scheme shown in Figure 4(c). With IK sub-carriers, 
IS fbis scheme requires a twofold increase for the correlation matrix size and FFT lengths, 
when calculating and y/j) respectively. 

Channel estimator conqjlexity can be further reduced (where the exponential 
power delay profile of the channel can be approximated as imiform), if the length of 
the uniform power delay profile is chosen correctly reduced complexity weighting 
20 coefBcients result. The length of the power delay profile is usually set to the cyclic 
prefix length. "Good'* CycKc Prefix (CP) length windows are (l+N/2) or (l+N/4), 
where N is the length of the IDFT used to form the OFDM symbol. In this way the 
real and imaginary parts of the fixed weighting matrix values are made to contain 
equal or zero entries when "good" cyclic prefix length windows are chosen. 

25 With a uniform power delay profile, coefBcients of the fixed weighting matrix 

can be made real if the Channel Impulse Response (CIR) is located symmetrically 
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around zero by perfoming a cycKc shift, as shown in Figure 5. This approach makes 
all the coefficients of the iSxed weighting matrix real, thus reducing flie complexity of 
the computations required to be performed by the channel estimator 28. 

Figure 5 (top) shows a typical channel impulse response 120. A 
5 uniform (rectangular) shaped power delay projBle 122 is drawn encompassing the 
impulse response. Figure 5 (bottom) shows both the chaimel impulse response and the 
assumed uniform power delay profile shifted to the left and therefore centering this 
power delay profile about zero. This is achieved by a cyclic shift when used with 
DFT/EDFT block processing, as used by OFDM systems. The negative time 
1 0 components appear at the end of the block as shown in Figure 5 (bottom). 

Returning once again to Figure 3, the sequence of steps carried out by tihie 
charmel estimator 28 in order to provide the LMMSE chaimel estimation by significant 
weigiht catching may optionally include the steps of performing, at step 92, a phase 
rotation of the least squares estimation matrix derived in stq> 82, and a complimentary 
15 step 94 of performing a de-rotation of the LMMSE channel estimates derived in step 
90. Finally, the channel estimation vectors are provided to the detector 26 in step 96* 

The cyclic shift for the chaime] impulse response can be achieved in the 
firequency domain by applying a linear phase rotation across the LS fi-equency 
estimates of (-27ikp/N), where the shift, p, is half ftie length of flie uniform power delay 
20 profile. Note p is negative for the complementary step of 94. Ihe latter step can be 
avoided if the data symbols are pre-rotated. 

If the "good" cyclic prefix windows are used, steps 92 and 94 may not be 
required. However, tiiis approach can reduce the results provided by the chaimel 
estimator 28 due to a less than optimal windowing of the channel impulse response. 

25 

The Applicants have carried out simulations in an 802.11a system with 2 
transmitters and 1 receiver. The mean squared error (MSE) for antenna J is given by: 

MSB, -^trac^£i^(Hj -Hj) (Hj ^jfj^ (12) 
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The system operated in an indoor HIPERLAN/2 non-sample-spaced channels 
A (r^ = 50 ns), B (r^ = 100 ns) aad C (r^ 150 ns), with the total transmit 
power normalized to unity. It was assumed that perfect knowledge of the SNR. and 
T^. were available for calculation of the nr,. 

5 The MSB channel estimation performance was evaluated by transmitting two 

long OFDM-BPSK pilot symbols through a fading multipafh channel 1000 times. For 
each iteration, the pilot symbols were simultaneously sent fix>m the two transmit 
antennas on interleaved sub-carriers. The duration of the two long pilots was 8 fis 
including double length CP of 1.6 yLs and the total system bandwidth was subdivided 
10 into IT = 52 sub-carriers (out of a possible 64). For the sparse approximations, the 
number of complex multipliers (M < K/2) was chosen to give targeted MSE error floor 
^25 dB. 

It was observed that the LMMSE by Single Value Deconq)osition (SVD) 
outperforms the LMMSE by Significant Weight Catching (SWC) in channel A, when 
15 the rank r ^8 as can be seen in Figure 6. At a fixed value of SNR = 25 dB, its MSE 
error floor is well below of 25 dB and the estimator requires 12 complex multipliers. 
However, if the charmers delay spread is increased (channels B and C), the LMMSE 
by SWC is a better compromise in performance versus complexity, as shown in Figure 
6. 

20 The LMMSE by SWC requires only 12 complex multipliers in order to reach 

an adequate performance in channel B and the estimator complexity is reduced by 
more than 50% compared to the fiill LMMSE. It should also be noted that the 
performance of the sinq>lified LMMSE algorithm remains almost unchanged in all the 
diannels, especially for the low number of complex multipliers (:^12). To illustrate 

25 the performance for a dynamic SNR range, the MSE in channel B is presented in 
Figure 7. The number of complex multipliers M - 3r/2 in the sparse s^roximations 
was set to the fixed nominal values of 12 and 21. With the MSE gain of 9 dB over the 
LMMSE by SVD for M = 12 at SNR = 3 0 dB, it can be seen that the LMMSE by SWC 
is the better choice for a reduced complexity LMMSE channel estimator. 

30 From the foregoing, it is apparent that LMMSE by SWC estimation technique 

described above can reduce computational complexity of the traditional LMMSE 
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chaimel estimator by more than 50% and it outperfbims the LMMSE by SVD when 
channel delay spreads exceeding 50 ns. 

Finally, it is to be understood that various modifications and/or additions may 
be made to the above described method of channel estimation vdthout departing from 
5 tibie ambit of the present invention as defined ia ihc claims appended hereto. 



